THE PRECIPITATION OF ENERGETIC HEAVY IONS INTO THE UPPER 

ATMOSPHERE OF JUPITER 

(KAS1-TB-101 155} TEE £ RECIPITATICN OF N88-24565 

E8EHGETIC HEAVY ICKS IBIO TEE UPPER 
ATHQSPHEEE OF JUPITER (BASA) 74 p CSCL 03B 

Unclas 

G3/9 1 0146661 

M. Horanyi* 

Supercomputer Computations Research Institute 
Florida State University 
Tallahassee, Florida 


T.E. Cravens 

Space Physics Research Laboratory 
Department of Atmospheric and Oceanic Science 
The University of Michigan 
Ann Arbor, Michigan 48109 


J.H. Waite, Jr. ^ 
Marshall Space Flight Center 
Huntsville, Alabama 


February 1987 


*On leave from Central Research Institute for Physics, Hungary 


Q3AI333H 



ABSTRACT 


Evidence for auroral particle precipitation at Jupiter was provided by the ultraviolet 
spectrometers onboard the Voyagers 1 and 2 spacecraft and by the International Ultraviolet 
Explorer (IUE). Magnetospheric measurements made by instruments onboard the Voyager 
spacecraft indicate that energetic sulfur and oxygen ions are precipitating into the upper 
atmosphere of Jupiter. We have constructed a theoretical model describing the interaction of 
precipitating oxygen with the Jovian atmosphere. The auroral energy is deposited in the 
atmosphere by means of ionization, excitation, and dissociation and heating of the atmospheric 
gas. Energetic ion and electron precipitation are shown to have similar effects on the 
atmosphere and ionosphere of Jupiter. 
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L INTRODUCTION 

The ultraviolet spectrometers on the Voyagers 1 and 2 spacecraft (Broadfoot et al M 
1981; Sandel et al., 1979) and the International Ultraviolet Explorer (IUE) satellite (Clarke et 
al., 1980; Yung et al., 1982) observed intense H 2 Lyman and Werner band emissions from the 

Jovian atmosphere at higher latitudes, and thus provided evidence for auroral particle 
precipitation at Jupiter. The Voyager ultraviolet spectrometer experiments (UVS) inferred 
about 60 kR of Lyman alpha emission and a total of about 80 kR of H 2 Lyman and Werner 

band emission if the latitudinal extent of the auroral zone was approximately 6000 km wide. 
This corresponds to an energy flux of 10 ergs/cm^/s in the auroral zone, or a total global power 
of 1-2 x 10*3 Watts if one identifies the precipitating particles as electrons (Waite et al., 1983; 
Atreya et al., 1981; Yung et al., 1982; Broadfoot et al, 1981; Gerard and Singh, 1982). Yung 
et al. (1982) interpreted IUE spectra of the Jovian aurora using a radiative transfer model for 
Lyman alpha and for the Lyman and Werner band emissions and they concluded that the 
precipitating electrons (if the particles are electrons) need to have energies between 1 and 30 
keV. Waite et al. (1983) presented a comprehensive theoretical model of both the auroral and 
non-auroral atmosphere and ionosphere of Jupiter and used this model to explore the 
aeronomical effects of auroral precipitation of 1 and 10 keV electrons. 

The Io plasma torus is thought to be the source of the energetic particles responsible for 
the observed aurora. Intense fluxes of charged particles have been detected in the inner 
magnetosphere of Jupiter by experiments on the Voyager spacecraft (Bridge et al., 1979; 
Krimigis et al., 1979; Vogt et al., 1979; Gehrels et al., 1981). Both the energetic and thermal 
ion populations in this region are largely composed of relatively heavy ions such as sulfur and 
oxygen. Most earlier magnetospheric studies of wave-particle interactions have suggested that 
electrons were more likely than ions to precipitate into the Jovian atmosphere (cf. Waite et al., 
1983). Goertz (1980) suggested proton precipitation and Thome (1981a, b) considered heavy 
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ion precipitation. 

Gehrels and Stone (1983) reported on observations of 1 to 20 Mev/nucleon oxygen, 
sodium, and sulfur ions in the Jovian magnetosphere by the cosmic ray subsystem (CRS) 
onboard Voyagers 1 and 2. Gehrels and Stone (1983) (designated as GS, hereafter) calculated 
the inward diffusion rate of these heavy energetic ions and assumed that they are scattered into 
the loss cone, thus resulting in auroral power input into the Jovian atmosphere. We reproduce 
in Figure 1 the auroral power spectrum shown in Figure 12 of GS. This figure shows the 
calculated total integrated power as a function of energy for oxygen and sulfur. GS also 
extrapolated these spectra to energies below those actually observed by the CRS in order to 
obtain the 1-2 x 10 13 W total auroral power required to explain the observed ultraviolet 
observations. 

Some additional evidence in favor of ion (rather than electron) precipitation comes from 
X-ray observations of Jupiter made by the Einstein X-ray observatory (Metzger et al., 1983). 
Metzger et al. (1983) showed that if Bremstrahlung radiation from electron scattering was 
assumed then a total auroral power of 10^ - 1()16 \y would be required to explain the 
observed X-ray emission. They suggested that K-shell excitation of energetic precipitating S 
and O could explain the X-ray observations with a reasonable total power in line with that 
required by the ultraviolet observations. 

This paper will investigate the nature of the energy deposition in the Jovian atmosphere 
by precipitating energetic oxygen ions and will explore the aeronomical effects associated with 
this deposition. The effects of energetic sulfur precipitation will not be considered explicitly; 
however, the general nature of the energy deposition should be independent of the exact atomic 
mass of a precipitating particle for comparable energies per nucleon. 

The questions that this paper will seek to address are: 

(1) How is the auroral power which is input at the top of the atmosphere distributed, or 
deposited, as a function of altitude? 

(2) How is this deposited energy divided among various aeronomical processes such 
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as ionization, dissociation, heating, and optical emissions, etc.? 

(3) How do these energy deposition efficiencies for various processes compare with 
the efficiencies calculated by Waite et al. (1983) for the precipitation of energetic electrons? In 
other words, how similar or different are the aeronomical effects of electron and heavy ion 
precipitation? 

(4) What are the optical emissions associated with the precipitation of oxygen and 
sulfur ions? Can these emissions be used for an observational test of the existence of heavy 
ion precipitation at Jupiter? 

Section 2 of the paper describes the physical processes included in both the energy 
deposition part of our theoretical model as well as the aeronomical part In section 3, we 
present the results of the energy deposition calculations. We present the results of the 
aeronomical calculations in section 4; in particular, we show calculated atomic hydrogen 
density profiles and calculated electron and ion density profiles. Section 5 is a summary and 
discussion section. And in a companion paper (Waite et al., 1987) we present the results of 
some IUE observations of the Jovian aurora. 
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2. BASIC PROCESSES 

The most important processes are charge-transfer, electron stripping, dissociative and 
non-dissociative ionization, excitation and elastic collisions. The competition between charge 
transfer and stripping determines the charge state of an oxygen ion at any given energy, but the 
bulk of the energy is lost via ionization processes. We have compiled the available data on the 
cross-sections for these processes, but in several cases we had to adopt various theoretical 
approximations to circumvent the lack of measured data. In our calculations we considered 
only four different charge states for the oxygen ions (0 £ q £ 3); the highest charge state is the 
triply ionized state. 

In this section we give a detailed description of the processes followed by a description 
of the construction of the stopping power function using these processes. The method of 
equilibrated beams is employed. Finally, we discuss the energy distribution of the secondary 
electrons produced by the stripping and ionization processes. 

2.1 The charge-transfer process 

The charge transfer process can be represented as: 

O q+ + X 0 (q ' 1)+ + X (2.1) 


where X stands for the target which can be H or H 2 . 

The cross sections for process (2.1) are plotted in Figures 2 and 3 for the different 
charge states of the oxygen ions in H and H 2 . The low energy cross sections are measured by 

Phaneuff et al. (1978), and the high energy cross sections by Phaneuff et al. (1982). For q = 1 
oxygen ions in H 2 the low energy part was measured by Nutt et al. (1979). The continuous 
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lines in Figures 2 and 3 follow the data and the dashed lines are "graphical extrapolations". 

Very little energy is lost by the oxygen during charge transfer collisions. Such 
collisions for q>l are exothermic, although we have assumed a small energy loss of 1 eV in 
order to account for a certain fraction of collisions possibly resulting in oxygen excitation in 
excess of the exothermicity. The energy loss for charge-transfer collisions of 0 + was taken to 
be the endothermic energy plus 1 eV (accounting for possible excess excitation of the product 
O atoms), giving a total energy loss of 2.8 eV. The energy available to the products of charge 
transfer collisions (equation 2.1) for q=3 and for a H 2 target is 40.5 eV (including the 

additional loss of 1 eV mentioned above) and most of this energy should be available to excite 
the product 0 + 2 ions so that O III emissions should be generated in a Jovian ion aurora. These 
emission features should appear below 1000A in the spectrum. Similarly, 20.7 eV is available 
for the q=2 charge-transfer process resulting in 0 II emissions. The exothermic energy for 
these collisions ultimately comes from the beam kinetic energy via electron stripping collisions, 
which are discussed in Section 2.2. 

2.2 The electron stripping process 

The stripping (or electron removal) process can be represented as: 

+ X -> O q+ + X + e (2.2) 


where X stands for H or H 2 . 

For this process we were unable to locate any measured cross section data for oxygen 
ions. The measured stripping cross-section of Goffe et al. (1979) for boron and carbon were 
very similar in both H and H 2 . We assumed that the oxygen cross sections are also similar, 

and we used the Goffe et al. (1979) data for the q=2 and for q=3 states. As there was no data 
for the q=l charge state, we constructed the cross-section for q=l using the theoretical work of 
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Karasima and Watanabe (1983). The resulting cross-section functions are plotted in Figure 4. 

The corresponding energy losses are the relevant ionization energies of oxygen (13.6 
eV for q=l, 35.1 eV for q=2 and 54.9 eV for q=3) plus the average energy of the secondary 
electrons (T). We assumed that the energy distribution function of the secondary electrons 
emerging from stripping (and in the reference frame moving with the oxygen ions) is similar to 
the distribution for ionization at any given energy of the incident oxygen ion. 

2.3 Ionization and excitation processes 

The ionization of H can be represented by: 

O q+ + H O q+ + H + + e (2.3) 


Shah and Gilbody (1981) measured the cross-section for this process for oxygen ions 
with 2 £ q £ 5 in the energy range of 12 £ E £ 195 keV/amu. Part of their data is plotted with 
thick lines in Figure 5. These measurements indicate that the ionization cross sections for 
oxygen ions with q > 1 are roughly independent of the charge state in the energy range below 
where the cross sections have their maxima; but at higher energies, where the cross sections 
have decreased below their maximum values, the cross sections scale with q as: 


To supplement these experimental results we used the calculations of Olson et al. 
(1978). They calculated the cross section for electron removal from H (that is, the sum of 
the cross sections of charge transfer and ionization, aj). They did this for ions with 1 < q 
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$ 50 in the energy range of 50 £ E £ 5000 keV/amu, and their results fall on a single curve 
given by the expression: 


a T = 4.6qA(l - exp (-A)) 10 16 


(2.5) 


where A = 32q/E, E has units of keV/amu and sj is in units of cm^. 

We assumed that the experimental trends evident for q > 1 were also valid for q = 1 
as well, and we used the measured data together with equation (2.5) to construct cross 
section functions for q = 1 and to extend the energy range for q = 2 and q = 3 cross section 
functions. 

We could not locate data in the literature for q = 0. Stripping and ionization are 
basically the same process if the projectile and target atoms are both neutral, and one calls it 
stripping or ionization depending on whether or not the target or the projectile is ionized. 
Hence, ionization and stripping cross sections for q = 0 should be very similar, we have 
assumed that these cross sections are the same. Information on the stripping of H in 
collisions with various heavier neutral targets was also used in trying to construct a realistic 
cross section (McNeal and Birely, 1973). The resulting cross section functions are plotted 
in Figure 5. The energy loss for these processes is the sum of the ionization energy (13.6 
eV) and the average energy of the secondary electrons (T). 

The ionization process for H 2 is more complicated than for H, as there are several 

possible final channels. For example: 


O q+ + H2 -> 

O q+ + h£ + e 

(2.6a) 


O q+ + H + + H + e 

(2.6b) 

O q+ + H2 -> 

O q+ + 2H + + 2e 

(2.6c) 
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The dominant channel is non-dissociative ionization. These ionization processes can be 
divided into two groups: the non-dissociative part (equation 2.6a) and the dissociative part, 
which is the sum of all channels leading to the dissociation of H 2 . The thick continuous 

dashed lines in Figure 6 represent the measured data for non-dissociative and dissociative 
ionization processes, respectively. The thin lines are extrapolations based on the same 
(perhaps questionable) arguments that we adopted for H. The q = 0 ionization cross 
section for H 2 was constructed in almost the same way as the H ionization cross section 

discussed above. 

The energy loss in the non-dissociative channel is the sum of the ionization energy 
(15.4 eV) and the average energy of the secondary electrons (T). The energy loss in the 
dissociative channel can be approximated as the sum of the dissociation energy of the H 2 

molecule (4.5 eV), the ionization of an H atom (13.6 eV) and the average energy of the 
secondary electrons (T). 

To estimate the cross section function for excitation of H and H 2 we assumed them 

to be proportional to the ionization cross-section for H and to the non-dissociative 
ionization cross section for H 2 , respectively. We assumed a proportionality factor of 0.6, 

independent of the incident energy and the charge state of the oxygen ion. The ratio of the 

total excitation cross section to the ionization cross section for electron impact on H 2 is 

approximately 0.6 at higher energies. The majority of this excitation is to the B^Z + and 

© 

C!n states resulting in Lyman and Werner band emission. 

2.4 The equilibrium fractions 

If the collision mean free path is less than the atmospheric scale height, then the 
population distribution of the different charge states in the incident ion beam will reach an 
equilibrium distribution at any given energy. Let us denote that fraction of oxygen ions 
with charge state q as <f>q, with the normalization condition: 
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5 >,=> (2 - 7) 

q 

Naturally, the <))q depend on the energy of the beam. The time evolution of the <J>q are 
governed by the equation: 


d<t> n V“' 

<2-8) 

q 

where Pqq 1 is the transition probability in unit time from the charge state q to q'. 

Assuming that charge distribution is in equilibrium (d<(>q/dt = 0) and that only processes 
with a change of one electron are important (Pqq* = P 8q', q+i), we find the conditions of 
the detailed balance: 


^q,q+l “ ^q+1 ^q+1, q (2*9) 

Equation (2.9) can be solved for each <j>q at any given energy by replacing the transition 
probabilities with the cross-sections for stripping (Pq,q+i) and charge transfer (Pq.q-i). 

The resulting equilibrium charge distributions are plotted in Figures 7 and 8 for H and H 2 , 

respectively. The results show that below approximately 10 keV/amu the majority of the 
oxygen particles in both H and H 2 will be neutral. With increasing energy, higher and 

higher charge states become dominant The highest charge state we considered was q = 3 
which actually includes all higher charge states as well (equation 2.3). Consequently, our 
calculated stopping power (discussed in subsection 2.6) becomes increasingly inaccurate 
for E>300 keV/amu. 


2.5 Elastic collisions 
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At high energies ions lose very little energy through elastic collisions with the nuclei 
of the medium, but this is no longer true at lower energies where the average net charge of 
the ion approaches zero. While the electronic part (inelastic collisions) of the stopping 
power is dropping to zero with decreasing ion energy, the average loss from the elastic 
collisions is increasing rapidly. 

We adopted the universal nuclear stopping-power function from Lindhard et al. (cf. 
Northcliffe and Schilling, 1970), which was based on the Thomas-Fermi interaction 
model. The nuclear and electronic (to be described later) parts of the stopping power are 
plotted in Figure 9. 

2.6 The stopping power 

We are primarily interested for the present study in the range of energy (1.0 £ E ^ 
300 keV/amu), which is well below the range of validity of the Bethe formula. The only 
information we were able to locate in the literature was the range and stopping power tables 
of Northcliffe and Schilling (1970) which does not contain any information about the 
individual processes and whose lowest energy point is at 12 keV/amu. Based on the 
individual cross-section data which we have compiled (and in some cases extrapolated) one 
can determine the stopping power function: 


s<E)= sr =s ‘ +s » (210) 

where n is the number density of the target medium. The stopping power S is the sum of 
the electronic stopping power (due to inelastic collisions) S e , and the nuclear stopping 
power (due to elastic collisions) S n . 


The electronic part of the stopping power is: 
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s . -XZfX AE * (2 - u) 

q » 

where and A Eqj are the cross section and energy loss, respectively, for oxygen ions of 
charge state q and for the i-th process. Replacing <J)q with their equilibrium values we can 
directly calculate S e> 

In Figure 9 we have plotted separately the two parts of the stopping power for H 2 

and, for comparison, the relevant stopping power from the Northcliffe and Schilling (1970) 
tables as well. The stopping power for H is very close to half of that for H 2 . Below 

approximately 200 keV/amu, our calculated stopping power is in reasonable agreement 
with the Northcliffe data, but at higher energies the difference is more significant. There 
are several possible reasons for the discrepancy. In the low energy regime we have more 
confidence in our results and despite the uncertainties of the cross section data we do not 
believe that the Northcliffe tables are any more reliable. At higher energies (E > 200 
keV/amu) , however, the buildup of charge states with q > 3 will introduce errors into our 
calculated stopping power, causing an underestimation of stripping and ionization cross 
sections which increases with charge state q. 


2.7 The energy distribution of secondary electrons 


To circumvent the lack of measured data for the energy distribution of secondary 
electrons produced in the ionization or stripping processes, we have assumed that the shape 
of the energy distribution function of the secondary electrons is independent of the charge 
state of the incident particle, and that it is the same for ionization and stripping in the frame 
of reference of the electron emitting species (i.e., H or H 2 for ionization and O for 

stripping). With these assumptions we can use the partial differential cross section of 



14 


secondary electrons for proton collisions with H 2 measured by Rudd (1979) to represent 
the shape of the secondary electron distribution for oxygen collisions on H and H 2 , These 
relative cross sections versus secondary electron energy were then integrated and 
normalized using the total ionization cross sections (integrated over all secondary electron 
energies), which were taken from Figures 5 and 6. 

We used the empirical fit of Rudd (1979) to the measured distribution function: 


o(Y,T) = 5rca 2 



NjlH CY7 IQ 
Ij(4 + (Y/I 1 ) 2 ) 


exp 



( 2 . 12 ) 


where T is the energy of the secondary electron in eV, Y = Ep/1836, Ep is the incident 
proton energy in eV, aQ is the radius of the First Bohr orbit (0.529 x 10'^), Ifj = 13.6 eV, a 
= 1.28 and Nj is the number of electrons in the target with binding energy Ij. We assume 
that Ep is the oxygen energy per amu. The electron production cross section for stripping, 
a (Y,T), is not given by equation (2.12), but is determined by a transformation of equation 
(2.12) to the oxygen frame of reference. This was done numerically. Figure 10 shows 
relative differential cross section for the stripping process as a function of electron energy 
for several primary energies. At low primary energies the relative cross section is almost 
the same as the ionization cross section (equation 2.12), but at higher energies the 
transformation to the oxygen reference frame shifts the cross sections to larger values of T. 

The agreement of equation (2.12) with the experimental ionization data is good 
below approximately 50 keV incident proton energy for all secondary electron energies. 

For proton energies above 50 keV, the low energy part of the secondary electron 
distribution is well represented by equation (2.12), but at higher electron energies the fit 
deviates from the data. For all incident oxygen energies we used equation (2.12) to 
calculate the average energy of the secondary electrons: 
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jTa(Y,T)dT 

T = -2- (2.13) 

J c(Y,T)dT 
o 

where I is the average ionization energy (= 15 eV) for H 2 . One finds that T = IY 1/2/2 for 
ionization, but is greater than this for stripping. The final value of T is determined by 
weighting the T for both ionization and stripping with the respective total cross sections. 

2.8 Basic Processes for Secondary Electrons 

The secondary electrons, produced by the interaction of the energetic ion beam with 
the H 2 and H in the upper atmosphere of Jupiter, will themselves interact with the neutral 

atmosphere. These electrons can ionize, dissociate, electronically excite, and vibrationally 
excite H 2 molecules, as well as heat the ambient thermal electrons in the ionosphere. Most 
secondary electrons for atmospheric densities greater than - 2 x 10^ cm* 3 (i.e., altitudes 
less than = 750 km) have mean free paths less than an atmospheric scale height and deposit 
their energy locally. But electron transport is important at higher altitudes. 

The two-stream method is used in our model to calculate secondary electron (and 
photoelectron) fluxes and the associated energy deposition rates for various processes. 

This method was described in Banks and Nagy (1970) and Nagy and Banks (1970). The 
pitch angle distribution of electrons is approximated by two streams, one going up parallel 
to the magnetic field, and the other going down. The energy bin width was 0.5 eV for 
electron energies less than 10 eV, 1 eV for energies greater than 10 eV and less than 100 
eV, and various larger bin widths were used for more energetic electrons. The electron 
impact cross section data for H 2 and H used in the model was discussed by Waite et al. 

(1983). 
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2.9 Aeronomical Processes 

Energy deposition drives a variety of aeronomical processes in the upper 
atmosphere of Jupiter. For example, dissociation of H 2 by energetic particles and their 

secondary electrons produces H atoms which diffuse down to the homopause region where 
they participate in hydrocarbon chemistry and eventually recombine. And ionization of H 2 

and H by energetic particles, or by solar radiation, leads to the formation of an ionosphere. 

For this paper, we used the aeronomical model described by Waite et al. (1983) and Waite 
(1981). 

The neutral temperature structure for the Jovian thermosphere was taken from 
interpretations of Voyager UVS stellar and solar occultation data by Atreya et al. (1981) 
and Festou et al. (1981). The neutral part of the model solves the continuity and 
momentum equations for H 2 , He, H, CH 4 , CH 3 , C 2 H 2 , C 2 H 4 , and C 2 H 5 , and also for 
the major ion in the ionosphere, H + . Photochemical solutions were obtained for C 2 H 5 +, 

CH 5 +, CH 4 +, CH 3 + CH 2 + , CH+, He+, HeH+, H 2 + , and H 3 + The ionospheric and 
neutral models were coupled with the energy deposition model so that consistent solutions 
are obtained. Ionization by solar extreme ultraviolet radiation was also included, as 
described in Waite et al. (1983); however, energy deposition by energetic particles (ions or 
electrons) is much more important than solar radiation in the Jovian auroral regions. 

The momentum equation for neutrals includes both molecular and eddy diffusion. 

We used an eddy diffusion coefficient of: 

K = 10 6 ( 10 13 /M) W cmV 1 (2.13) 

where M is the total atmospheric number density. This choice of K gives reasonable 
agreement between calculated densities of CH 4 , C 2 H 2 , and C 2 H 6 near the homopause and 

hydrocarbon mixing ratios inferred from Voyager UVS stellar occultation results (Atreya et 



17 


al., 1981) 

Both the neutral temperature structure and the eddy diffusion coefficient were 
derived for low latitudes, and not for the higher latitude auroral regions. Auroral heating 
and the associated thermospheric circulation will certainly alter the neutral atmosphere at 
higher latitudes. However, to a large extent the energy deposition (and even the associated 
aeronomy) of energetic precipitating particles does not greatly depend on the details of the 
atmosphere structure if one uses atmospheric density level, rather than altitude, as the 
independent variable. In this paper, we will often include in the relevant figures both an 
altitude scale and a scale for total neutral number density (which is almost equal to the 
molecular hydrogen density except at the highest altitudes where H is relatively abundant). 
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3. ENERGY DEPOSITION 


The total power delivered to the Jovian atmosphere by energetic oxygen and 
sulfur ions as a function of energy was determined by Gehrels and Stone (1983)(see Figure 
1). We assumed that the power for both oxygen and sulfur was all in the oxygen ions 
since our cross section data is for oxygen. We also assumed that the energy spectrum 
could be fit with a power law which was also applied to energies much less than those 
which were actually measured (Gehrels and Stone, 1983). When the total energy flux is 
transformed to a particle flux (per unit energy) one obtains a E'3-l power law where E is 
the energy. Fluxes per unit area were obtained by assuming that the precipitating ions were 
incident on auroral zones 6000 km wide, centered at 65° latitude in both hemispheres. 48 
discrete beams with incident energy, E^ (k=l,48), were used to represent the incident ion 

energy spectrum from 0.1 keV/amu up to 1 Mev/amu (see Table 1). The energy bin width 
for beam k is AE^ = (E^+i - E^. i)/2. The net downward incident particle flux for beam k 

is: 

F k = 2 x 10 4 E k 3-1 A E k (cm* 2 s' 1 ) (3.1) 

where E^ and AE^ are in units of MeV/amu. 

The calculations in this paper used equation (3.1) with a range of lower energy 
cutoffs, Emin. However, the paper emphasizes the calculation for Ejjun = 40 keV/amu (or 
k>15); that is, we assumed zero input flux for E less than E^n and used equation (3.1) for 
higher energies. The total power per unit area for this choice of E^n is 15 ergs/cm^/s 
which corresponds to a total global power of approximately 2 x 10^ W. E^n = 40 
keV/amu was chosen so that the calculated intensity of the H 2 Lyman and Werner band 
emissions was approximately the 80 kR or so observed by the UVS experiment on VI and 
V2 (Broadfoot et al„ 1979, 1981). The total power per unit area required for precipitating 
10 keV electrons to produce the same H 2 Lyman and Werner band emission is 10 
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ergs/cm^/s (Waite et al., 1983). 

The downward pitch angle distribution was assumed to be isotropic and we put 
all the flux in at an "average" pitch angle with <cosa> = 0.5. The energy of an ion (or an 
ion beam k) as a function of altitude can be found from its initial energy at the top of the 
atmosphere by using the stopping power 


z 

E(z) = E 0 -2 J S [ E (z') ] n (z') dz’ (3.2) 

®p 

where z t0 p = 3350 km, n(z) is the neutral density at z, the stopping power S is from 
equations (2.9) and (2.10) The energy at the top of the atmosphere, E 0 , is equal to E^. 

Actually, we included the densities and stopping power for both H 2 and H, although only 
H 2 is important except at higher altitudes. The factor of 2 in equation (3.2) accounts for 
ratio of pathlength increment to altitude increment for the average pitch angle used. Note 
that the energy refers to the energy of an equilibrated beam composed of several charge 
states. 

The total energy deposition rate at a given altitude (in units of eV/cm^/s) for 
beam k, Qk(z), can be found from the energy for beam k at altitude z, E(z), and from the 

stopping power: 

Q k (z) = 2 F k n(z) S[E(z)] (3.3) 

where the factor of 2 is from the averge pitch angle. Again, both H 2 and H were included. 

The energy deposition rate at z for a complete spectrum is found by adding up the 
appropriate Qk(z) for the desired range of k. The energy deposition rate for a particular 

process, rather than the total deposition rate, is found by using only the relevant part of the 
stopping power S in equation (3.3). 
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3.1 Total Energy Deposition Rates Versus Cutoff Energy 

Total energy deposition rate profiles are shown in Figure 9 as functions of 
altitude (an atmospheric density scale is also provided) for a range of cutoff energies, 

Emin- The neutral atmosphere model used for these calculations is called "unconverged" 

which means that the H densities are small and only take into account solar processes 
(Waite et al., 1983). This atmosphere is mostly H 2 up to an altitude of ® 900 km. The H 
densities for the "converged" atmosphere are much larger and include in a self-consistent 
manner the H produced by the energetic ion precipitation. The maximum energy deposition 
rates in Figure 9 increase with decreasing values of E^n since more energy is input at the 

top of the atmosphere. And the altitude of this maximum deposition also increases with 
decreasing values of E ^in since less energetic oxygen does not penetrate as deeply into 

the atmosphere. The lowest energy in the Voyager energetic ion data is somewhat above 
300 keV/amu (Gehrels and Stone, 1983). The deposition rate profile for this cutoff has a 
peak at an altitude of 400 km (total number density of = 10 13 cm* 3 ) and has a total 
integrated power which is only about 10% of that required to explain the observed UV 
emissions from H 2 . The homopause and the top of the hydrocarbon layer are both located 
approximately at 350 km in the equatorial region, but are possibly located at higher 
altitudes for higher latitudes, so that ions more energetic than 300 keV/amu will penetrate 
into this layer. 

3.2 Energy Deposition Rate For 40 keV/amu Cutoff Energy 

A total energy input of 15-20 ergs/cm 2 /s is required to generate the observed UV 
airglow intensity, and in order to obtain this input from the power spectrum, a cutoff 
energy, E mj n = 40 keV/amu is required. This requires an extrapolation in the observed 
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energetic ion spectrum by about a factor of 10 in both particle energy and total power 
(Gehrels and Stone, 1983). Figure 10 shows the energy deposition rate profile for this 
cutoff energy. The peak deposition rate is 1.1 x 10 6 eV/cm 3 /s at an altitude of 430 km 
(n tot =3.8 x 10 12 cm* 3 ). The total integrated energy deposition rate is 15 ergs/cm 2 /s, or 
9.25 x 10 13 eV/cm 2 /s, giving a global power of 2 x 10 13 W. 

Some breakdown of how the beam energy is deposited in the atmosphere is also 
indicated in Figure 12. The ionization and the charge-transfer/stripping processes account 
for most of the deposited energy. About 33% of this energy is put into secondary 
electrons; and the rest goes into ionization of molecular hydrogen and excitation of the 
energetic oxygen. Excitation of H 2 is not a major loss process for the oxygen beam but 
still accounts for about half of the required UV emissions. We assume that the vast 
majority of H 2 excitation leads to Lyman and Werner band emission either directly or via 
cascade from higher-lying allowed levels. Nuclear/elastic collisions are not a very 
important energy loss process, and this energy ends up as part of the total neutral heating 
rate as will be discussed later. 

The fractional contributions of each charge state in the equilibrated beam to the 
energy deposition are shown in Figure 13 as functions of altitude. Energy deposition by 
0 + (q=l) and 0‘ H tq=2) dominates near the top of the atmosphere, as might be expected 
from the equilibrium fractions for energies in the 40 -100 keV/amu range (Figurs 7 and 8). 
As the precipitating ions penetrate deeper into the atmosphere and lose energy, the 
equilibrium fractions of q=0 and q=l become larger and these charge states dominate the 
energy deposition. However, very deep in the atmosphere the contributions of the q=2 and 
3 states begin to pick up again because the lower energy part of the input spectrum has 
completely degraded and only more energetic ions remain. 


3.3 Energy Deposition by Secondary Electrons 
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The secondary electron distribution (see section 2) is strongly weighted towards 
lower energies. Figure 12 shows the production rate of electrons as a function of electron 
energy, T, for an altitude of 450 km and for the input spectrum with Eujj n = 40 keV/amu. 

The shape of the production rate curve at other altitudes is very similar to this one, although 
the total altitude-integrated electron production rate is just equal, at each altitude, to the total 
electron production rate. The total electron production from the primary beam is due both 
to direct ionization and to stripping collisions as shown in Figure 12. Total H 2 + and H + 
production rates (including ionization by both the primary beam and by secondary 
electrons, and including the contribution due to the charge-transfer process) are shown in 
Figures 15 and 16. 

The two-stream method was used to calculate the upward and downward 
electron fluxes (0 + and <{>-, respectively) as a function of electron energy and altitude. 

Figure 15 displays <{> + and <}>* versus altitude for a typical electron energy of T = 9.25 eV. 
Notice that <}> + = 9* except for altitudes greater than * 750 km, well above the altitude of 
peak energy deposition. This is an indication that the electrons degrade and deposit their 
energy in the near vicinity of where they were produced and that electron transport is not 
important Some electrons do escape at higher altitudes where the mean free path becomes 
long compared to an atmospheric scale height. Figure 18 shows the differential electron 
flux as a function of energy for an altitude near the peak of the energy deposition (z = 450 
km). The shape of this distribution reflects both the shape of the electron production 
energy distribution (i.e., Figure 14) and the energy dependence of the electron impact 
crosssections. 

Electrons with energies greater than the ionization potentials of atomic or 
molecular hydrogen can generate further secondary ionization. Figures 15 and 16 also 
show the contribution of secondary ionization processes to the total ionization rates. 
Secondary electrons also electronically excite H 2 and H (producing Lyman and Werner 
band emission and Lyman alpha emission), and vibrationally excite and dissociate H 2 . 



ORIGINAL PAGE IS 
OE POOR QUALITY 

The altitude dependence of these various processes is very similar to that of seconds 
ionization (Figure 15) and will not be shown. The altitude integrated rates (that is, c 
production rates) for all these processes will be discussed in the next section. 

3.4 Column Production Rates and Energy Deposition 

We have integrated the production rates over altitude for all important pr 
for oxygen precipitation with Emin = 40 keV/amu. The contributions of the primar 

beam and the secondary electrons are distinguished. The fraction of the total input < 
deposited in the atmosphere via a particular process - as determined by the column 
efficiency -- is also calculated. Table 2a lists these results for the unconverged atm< 
(low H densities) and Table 2b for the converged atmosphere (large H densities). 1 
energy input is 15 ergs/cm 2 /s for the 40 keV/amu cutoff energy and the total calcula 
Lyman and Wemer band intensity is 66 kR which is in rough agreement with the V 
UVS results. The contributions of the primaries and the secondaries are about the 
Lyman and Wemer band excitation. The H 2 band excitation for the converged atn 
is about 20% smaller than for the unconverged atmosphere because more energy i 
deposited in H for the former than for the latter. On the other hand, the converge< 
atmosphere has much more Lyman alpha emission than does the unconverged ao 
Ionization is the most important process, accounting for about 45% of 
energy input. "Ionization energy" here refers to the production rate times the ior 
potential and does not include the energies of the secondary, or teritary, electror 
accounted for seperately. H 2 + production is larger than H + production for both 
atmospheres. Most of the energy put into H 2 + production eventually ends up a 
heat due to ion-neutral chemistry and ion recombination (Waite et al., 1983). r 
beam is much more important than secondary electrons for the ionization proc 
Secondary electrons are more important than the primary ion be air 
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processes such as direct dissociation of H 2 (resulting in H production), direct vibrational 
excitation of H 2 , and electron heating. For these processes, electrons with energies below 
=■ 12 eV make the greatest contribution. However, vibrational excitation can also result 
from the cascading associated with Lyman and Werner band emission. Similarly, some H 
production and some neutral heating is associated with H 2 + production because this ion 
typically reacts with H 2 in the Jovian atmosphere to produce H 3 + and H plus some excess 
energy. Most of the H 3 + then recombines dissociatively providing H and heat. Similarly, 
almost all of the energy put into vibrational excitation, electron heating, and dissociation 
also ends up eventually as neutral heat due to various quenching, recombination, and 
cooling processes. The net total neutral heating efficiency is close to 47% for the 
unconverged atmosphere, and is a couple percent less than this for the converged 
atmosphere. 

A significant pan of the energy deposited by the energetic oxygen via the 
electron stripping process reappears during charge-transfer collisions as excitation energy 
of the product energetic oxygen ions or atoms (10%). We estimate that * 5% of the total 
input energy ends up as excitation of oxygen in its various charge states. (The other 5% is 
attributed to sulfur.) This excitation should appear as optical emissions mainly from 01, 
OH, and Om with a total intensity over all spectral features of roughly 40 kR. Most of 
these spectral features are below 1000 A in the spectrum and should include spectral 
features observed by the Voyager UVS from the Io plasma torus. But since we have 
assumed that energetic sulfur acts like oxygen in its energy deposition and have only 
performed the calculations for the latter, there should really only be about 20 kR of excited 
oxygen emission and the remaining 20 kR or so should come from excited sulfur species. 
Furthermore, some incident oxygen and sulfur are also probably excited directly via 
collisions with the ambient H 2 ; however, we have no information concerning this type of 
process and have neglected it, assuming it is less important than excitation via the 
stripping/charge transfer process. This is discussed further later in the paper. 
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The miscellaneous category in the Tables represents: (1) energy put into 
processes not explicitly accounted for, such as rotational excitation, (2) errors in assigning 
specific energy losses for the various processes, (3) and errors in the calculation. For 
example, there is a several percent error due to not accurately keeping track of the beam 
energy in the last (and lowest) altitude bin for a given beam. Kozyra et al. (1981) showed 
that almost all of this energy ends up as neutral heat, so that one can consider that a large 
part of the energy in the miscellaneous category is neutral heat 

One of the goals of this paper is to compare the atmospheric effects of electron 
and ion precipitation. Table 3 is a comparison of the column energy efficiencies for ion 
precipitation (Table 2a) and 10 keV electron precipitation (from Table 5a of Waite et al., 
1983). The peak of the energy deposition profile for the 10 keV electron aurora is at 460 
km, almost the same altitude as the peak of the deposition rate for the oxygen beam with 
Emin = 40 keV/amu. The efficiencies for H 2 + production and for assorted "low energy" 

processes (i.e., direct vibrational excitation, electron heating, direct neutral heating, direct 
H production) are very similar for the two types of auroral precipitation. However, oxygen 
precipitation is more efficient than electron precipitation for either H + production or oxygen 
(and sulfur) excitation. But electrons are more efficient than the oxygen beam for Lyman 
alpha excitation and Lyman and Werner band excitation, which is why an energy input of 
15 ergs/cm 2 /s, or more, energy input is required for the oxygen beam and only 10 
ergs/cm^/s for electron precipitation. However, the error in the calculation of the 
emission rates for the oxygen beam is roughly equal to this difference. Overall, the energy 
deposition patterns for the two types of auroral precipitation are rather similar, especially 
taking into account the uncertainties in both sets of calculations. 
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4. AERONOMICAL EFFECTS OF ENERGETIC ION PRECIPITATION 

It was established in the last section that the pattern of energy deposition for electron and 
heavy ion precipitation is about the same, within a factor of two or so. The aeronomy 
associated with auroral precipitation depends on where in the atmosphere the energy is 
deposited as well as on the total amount of energy deposited. The altitude of the peak 
energy deposition for the precipitation of magnetospheric oxygen with energies in excess of 
40 keV/amu is about the same as for the precipitation of 10 keV elections (z = 460 km 
which is equivalent to a total atmospheric density of 3 x 10*2 cm*3), and, consequently, 
the aeronomical effects of the two types of auroral precipitation are similar. Waite et al. 
(1983) demonstrated that both the neutral atmosphere and the ionosphere are greatly altered 
by electron precipitation, and in this section we will show that this is also true for oxygen 
precipitation. 

4.1 Neutral Atmosphere 

The total column heating efficiency for the neutral thermospheric gas should take into 
account: (1) the direct heating due to the dissociation of molecular hydrogen by electron or 
ion impact, (2) electron heating because this heat ultimately goes to the thermospheric 
neutrals because the electrons cool via electron-neutral collisions, (3) vibrational excitation 
because virtually all of this energy ends up as thermal energy due to the collisional de- 
activation of the vibrationally excited H 2 , (4) chemical heating associated with chemical 
reactions initiated by H 2 + . Recombination of H can be neglected as a heat source above 
the homopause. Waite et al. (1983) found that the total neutral heating efficicency for 
electron precipitation was very close to 50% for both 1 keV and 10 keY electrons. We 
calculated this total heating efficiency for ion precipitation by adding up the relevant entries 
in Table 2, and we obtained a value of 47%, which is the same as the electron precipitation 
value within the error. 

The neutral thermospheric temperature profile can be determined by solving the energy 
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equation including heat conduction, heating processes, and cooling processes. The neutral 
temperature calculated in this manner by Waite et al. (1983) was 1000-2000 K higher for 
the auroral case than for the case with solar EUV heating only. We will not repeat this 
calculation for the ion aurora because the resulting temperature profile would be the same 
for equal auroral energy inputs, because the heating efficiency is the same for electron and 
ion precipitation and because the altitude of the maximum heating is almost the same (for 10 
kev electrons). The temperatures one would calculate for an ion aurora would be 
somewhat larger than for an electron aurora, because the Lyman and Werner band 
excitation efficiency is about a factor of two less for ions than for electrons, so that a larger 
incident power is required to obtain the observed intensities for the former. 

Neutral composition is also affected by auroral energy deposition. In particular, the 
abundance of atomic hydrogen should be very large in auroral regions because dissociation 
of molecular hydrogen, both directly by particle impact and indirectly due to photochemical 
reactions initiated by H2 + , is greatly enhanced in auroral situations. Altitude profiles of H 2 

, H, and He densities are displayed in Figure 19. The curve labeled "unconverged" 
designates the H density profile calculated without any auroral precipitation; the peak 
density is about lO 1 ^ cm'3. The H profile labeled "converged" was calculated including 
the auroral production of H. The composition and the precipitation/energy deposition parts 
of the model were iterated until convergence was achieved. The peak H density for the 
auroral case is 4 x 10^ * cm*3; this is very close to the peak density obtained by Waite et al. 
for 10 keV electron precipitation. Both sets of auroral calculations ignore the possible 
dilution of H in the auroral regions due to horizontal transport of H out of these regions via 
the thermospheric winds generated by the auroral thermospheric heating. 

Vibrationally excited molecular hydrogen is another component of the neutral 
atmosphere strongly affected by auroral precipitation (Cravens, 1987). But again, the total 
energy input and location in the atmosphere is more important than the species of the 
precipitating panicles. The column energy efficiency for vibrational excitation is about 
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20% for either incident electrons or ions. Cravens (1987) calculated the distribution of 
vibrationally excited for several scenarios including the 10 keV auroral precipitation case 
considered by Waite et al. (1983). It is reasonable to assume that if these vibrational 
calculations were re-run using ion aurora inputs, then the results would be similar to the 10 
keV electron aurora results. 

4.2 The Ionosphere 


Ionization rates are greatly enhanced in the auroral regions and it was shown in the last 
section that this is true for ion as well as electrons. The production efficiency for H 2 + is 

almost the same for the two types of aurora, although the peak production rate for ion 
aurora must be somewhat larger due to the somewhat larger incident power needed to 
produce the observed intensities of Lyman and Wemer band emissions. However, the 
efficiency for H + production is 3 times greater for ion precipitation than for electrons. In 
both cases, H + is the major ion even though the production rate of H2 + is much higher 

than that of H + . The reason is that for standard models in which radiative recombination is 
the most important loss process for H + , the chemical lifetime of this ion is much longer 
than the lifetime of other ions such as H 2 + and H 3 + (c.f., Atreya and Donahue, 1976). 

H 2 + is rapidly converted to H 3 + via the following reaction: 

H/+H 2 ^H 3 + +H (4.1) 

Then the H 3 + ions are rapidly removed via dissociative recombination: 


H + + e -> H + H 

3 2 


(4.2) 


-» H+ H + H 
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The density of H 3 + depends both on the production rate of H2 + and on the rate coefficient 
for reaction (4.2); we have used the rather s mall value of 10*8 cm^ s'* in agreement with 
recent theoretical and experimental results (Michels and Hobbs, 1984; Smith and Adams, 
1984; MacDonald et al., 1984). 

The calculated density profiles for H 3 + and H + , as well as electrons, are shown in 
Figure 20. The densities of H 2 + and other molecular ions are extremely small and are not 
shown. H 3 + becomes the major ion only at altitudes well below the ionospheric peak. H + 
is the major ion at higher altitudes and the maximum electron density is 3 x 10*7 cm*3. The 
peak electron density for thelO keV electron aurora was only 1.5 x 10? cm* 3 (Waite et al., 
1983), which is reasonable considering that: (1) the peak H + production rate for the ion 
aurora is about 4 times larger than it was for the electron aurora and, (2) the electron 
density varies as the square root of the production rate for a major ion in photochemical 
equilibrium. In any case, the observed electron densities at higher latitudes on Jupiter 
(Eshleman et al., 1979) are about a factor of 100 smaller than the calculated values. 

Clearly, some additional loss mechanism is required for the H + ions and several 
mechanisms have been suggested (cf., Atreya et al, 1979; Waite and Cravens, 1986) 
including the reaction with vibrationally excited H 2 with v £ 4: 

H + + Hj (vS4) -» h/ +H (4.3) 

Cravens (1987) demonstrated that this can be a potent loss process in the auroral region. 
This result should be as applicable to an ion aurora as to an electron aurora. 

Connemey and Waite (1984) suggested another loss process that might be operating as 
a sink of H + ions on Saturn. Water vapor entering the top of the atmosphere from the 
rings could act as a chemical sink, converting ionospheric protons to H 20 + ions which 

have a short lifetime. A similar process might operate in the Jovian ionosphere if heavy ion 
precipitation is occuring. The energetic O ions/atoms eventually neutralize and thermalize. 
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They will then diffuse downward. If the thermal oxygen density is sufficiendy large then 
H + ions could be removed via a charge exchange reaction: 

H + + O -> 0+ + H (4.4) 

The reaction rate coefficient for this is about 10'9 cm3 S -1 (Banks and Kockarts, 1973) and 
an O density of roughly 10^ cm' 3 near the ionospheric peak would be sufficient for this 
reaction to compete with radiative recombination as a removal mechanism for H + . The 0 + 
ions are rapidly removed by reaction with H 2 molecules: 

0 + + H2 -> OH + + H (4.5) 

OH + will recombine dissociatively reforming the original 0 atom We are currently 
looking into this mechanism in more detail. 
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5. K-SHELL X-RAY EMISSION FROM OXYGEN AND SULFUR 

X-ray observations of the Jovian aurora were made by the Einstein X-ray Observatory 
(Metzger et al., 1983). Metzger et al. (1983) estimated that a total auroral power of 10^ - 
1016 w would be required to explain the observed emissions if they were due to 
Bremsstrahlung associated with energetic electron precipitation. However, they estimated 
that a total power of only 10^ - 10^ W would be required to explain the observations 
with K shell excitation of the precipitating oxygen and sulfur ions. The X-ray photons for 
K shell excitation of oxygen and sulfur are 540 eV and 2300 eV, respectively. The total 
power required to explain the observations in this way was 5 x 10^ W at the lower energy 
and 4-8 x 10^ W at the higher energy (Metzger et al., 1984). 

We have calculated K shell emission rates using our model. Our model is only 
applicable to oxygen but we can assume that the pattern of energy deposition for 
precipitating oxygen and sulfur is very similar and apply the model to both species. We 
used the K shell excitation cross sections of Garcia et al. (1973) and for the x-ray emission 
efficiencies of sulfur and oxygen, we used 10% and 0.8% respectively. The K-shell 
excitation cross sections are very small for ion energies less than 100 keV/amu, and they 
are proportional to q^ where q is the charge on the ion. As a result, only the more 
energetic part of the incident auroral beam contributes to K-shell excitation. 

The calculated energy efficiencies (relative to the total auroral power for both oxygen 
and sulfur, and assuming equal energy inputs for both species) for X-ray emission are 4 x 
10'6 and 8 x 10‘7 for oxygen and sulfur, respectively. Using these efficiencies, one finds 
that the total power required to explain the X-ray observations is 0.7 - 1.5 x 10^ W. 
However, there is a serious source of error in these figures. Our model is quite inaccurate 
for energies in excess of 200 keV/amu or so because we terminate our equilibrium charge 
state at q=3, although at higher energies the actual values of q can be significantly larger. 
This is not a serious problem for the overall energy deposition since 90% of the incident 
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energy flux is for ions with energies less than approximately 200 keV/amu; however, it is a 
serious problem for K-shell excitation. We estimate that the "average" value of q at higher 
energies is 4-6, and the corrected K-shell excitation cross section then becomes about 4 
times larger than we used. The revised total power required to explain the X-rays is now 
about 2 x 10 14 W. This value refers to all incident ions in excess of 40 keV/amu. 
However, if one were to truncate the energy spectrum at about 200 keV/amu, then the X- 
ray excitation rate would be affected very little and the total auroral power required would 
be 1*3 x lO 1 ^ w. This auroral power is a factor of ten greater than that deduced by 
Gehrels and Stone (1983) from the magnetospheric measurements for these incident ion 
energies; however, there are uncertainties in the deduced auroral energy input, in the X-ray 
observations, and in our calculations of the K shell excitation efficiencies. 
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6. ULTRAVIOLET EMISSIONS FROM PRECIPITATING OXYGEN AND SULFUR 

The oxygen ions/atoms in the incoming auroral beam can be directly excited via 
collisions with the atmospheric gas. 

O q + Hj O q+ * + ^ (6.1) 

OQ+ can also be excited during a charge transfer collision: 

0 (q+I)+ + Hj O q+ * + h/ (6.2) 

Reaction (6.2) is highly exothermic for q > 1, and the product oxygen ion is likely to be 
highly excited and/or the H2 + ion dissociated. See the companion paper by Waite et al. 

(1987) for additional discussion of these processes. These excitation mechanisms also 
apply to sulfur. 

Two emission features in particular which should be present in the auroral spectra 
measured at Jupiter are the 01 1304 A line and the SII 1256 A line. Waite et al. (1987) 
reported on IUE observations of the Jovian aurora in the 1240 - 1500 A spectral region. 

The measured intensities at 1256 and 1304 A were at the noise level, and upper limits of 
2kR and 150 R, respectively, were set for these two features. First, consider the direct 
excitation mechanism. Waite et al. estimated that the excitation cross section of O colliding 
with H 2 was 4% of the electron removal/stripping cross section for O. Our calculations 

show that the column rate of electron stripping for O is 5 x 10^ cm'2 s"^ if the entire 
incident beam is composed of oxygen and 2.5 x 10^ cm*2 s" 1 if the beam is half oxygen 
and half sulfur. The estimated OI 1304 A intensity is 1 kR if one multiplies this latter 
number by the 4% emission efficiency. 1 kR is a factor of 10 larger than the upper limit set 
by the IUE measurements. 

Now consider the charge-transfer excitation mechanism. The total column energy 
deposition efficiency (for all charge states) that is available for exciting the product 
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ions/atoms from the charge transfer process is close to 10% (see section 4). The 
corresponding column energy dissipation for this process is =10^ eV/cm^/s for a total 
auroral power input of =* 10*3 w. Following the discussion given by Waite et aL (1987) 
we can make the following estimates of O and S spectral intensities. Suppose that half of 
the available energy goes into dissociating the product H 2 + ions and half into exciting the 

product O and S ions/atoms. The total column emission rate for all S and O features is then 
= 20 kR for a typical UV photon energy of 20 eV or so. There are approximately 10 strong 
emission features in the UV spectra of these two species, so that one can expect that 
roughly 2 kR is emitted per spectral feature (on the average) for SI, SII, SHI, SIV, and 
OI, OH, and Om (Waite et al., 1987). None of these features in the IUE spectral range 
(and in particular, SII 1256 A) have intensities much above the noise level in the IUE 
spectra. 
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7. DISCUSSION 


One goal of this paper, and of the companion paper by Waite et al., was to ascertain 
whether or not auroral ion precipitation is occurring on Jupiter. Another goal was to 
determine how heavy ion precipitation affects the upper atmosphere and ionosphere of 
Jupiter. Our basic conclusion concerning this second point is that the energy deposition 
and the associated aeronomy are very much alike for both electron and ion precipitation. 

For example, about 30-40% of the incident auroral power goes into ionizing H 2 for both 

types of auroras. However, there are also some differences. Energetic ions are three times 
more effective than electrons in producing H + ions (if one extrapolates the Gehrels and 
Stone spectra to lower energies). And energetic electrons are approximately twice as 
effective as ions (14% versus 7%) for the excitation of the H 2 Lyman and Wemer bands. 

One consequence of this latter difference is that a column energy input of 20 ergs/cm^/s is 
needed for ion precipitation to generate 80 kR of these features, whereas only 10 
ergs/cm^/s is needed for electrons. The input for our model was only 15 ergs/cm^/s, and 
thus only 66 kR was produced. These are minor differences, though; basically, the two 
types of precipitation have similar effects on the atmosphere. The H densities and the 
electron densities calculated for the two types of particle precipitation are the same within a 
factor of 2. 

We will now briefly discuss how likely ion precipitation is at Jupiter. A more complete 
discussion of this issue is provided in the companion paper. Measurements of energetic 
sulfur and oxygen in the Jovian magnetosphere which were made by instruments on the 
Voyager spacecraft (Gehrels and Stone, 1983) suggested that these ions were being lost via 
pitch-angle scattering and subsequent auroral precipitation of these ions. The X-ray 
observations made by the Einstein Observatory also support the existence of ion 
precipitation (Metzger et al., 1983). However, other magnetospheric evidence suggests 
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that electron precipitation is also taking place (cf., Waite et al. t 1987). The lack of 
observed O and S emissions in the UV also throws some doubt on the degree of S and O 
precipitation. 

The following scenario might be able to reconcile the conflicting evidence concerning 
ion versus electron precipitation. Magnetospheric O and S does precipitate, but only for 
energies in excess of 100 - 300 keV/amu; that is E m ^ n = 200 keV/amu. This does not 

contradict the Voyager measurements of Gehrels and Stone (1983) since the ions they 
observed exceeded even this energy. The total input power of energetic ions is now about 
2 x 10*2 w and the beam will penetrate deep into the atmosphere (see Figure 11). The 
excitation rates for the H 2 Lyman and Werner bands and for the S and O UV features are 

now all a factor of 10 smaller. The UV intensities are even smaller if the beam penetrates 
the hydrocarbon layer where UV photons can be absorbed. The predicted S and O line 
intensities will now be in line with the upper limits set by the IUE observations. Of course, 
the Lyman and Werner band intensities will also be much reduced, and will now be less 
than 10 kR; some other mechanism (such as electron precipitation) is now required to 
produced the observed intensity of = 80 kR. The K-shell emission rates for O and S 
remain almost unaffected by the increase in E min ; they are still somewhat too small, but 

one is no worse off with the new scenario than one was with the "full" ion aurora. A 
combination of uncertainties might be able to explain the discrepancy between the model 
and observed X-ray intensities. 

In summary, we are suggesting that: (1) energetic S and O ions with energies in excess 
of 100-300 keV/amu, and with a total power of * 10*2 w, are precipitating into the upper 
atmosphere of Jupiter and that these ions can partially or entirely explain the observed X- 
ray emissions, and (2) energetic electrons with incident energies in excess of 10 keV and a 
total power of = 10^ W are also precipitating and explain the bulk of the observed Lyman 
and Werner band emissions. Atmospheric and ionospheric processes are not very 
sensitive to the identity of the precipitating particles (i.e., electrons, heavy ions, or protons) 
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but depend only on the total amount of energy deposited in the atmosphere and on the 
atmospheric level at which this energy is deposited. 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


Figure 7 


Figure 8 


Total integrated power for precipitating oxygen and sulfur ions versus 
oxygen/sulfur energy. From Gehrels and Stone (1983). 

Charge transfer cross sections versus oxygen energy for oxygen ions of 
different charge states in atomic hydrogen. 

Charge transfer cross sections versus oxygen energy for oxygen ions of 
different charge states in molecular hydrogen. 

Electron stripping cross sections versus energy for oxygen ions of different 
charge states in atomic and molecular hydrogen. 

Cross sections for the ionization of atomic hydrogen versus energy for 
different oxygen charge states. 

Cross sections for the ionization of molecular hydrogen (dissociative and 
non-dissociative) versus energy for different oxygen charge states. 

Equilibrium fractions of an oxygen beam in various charge states versus 
energy for atomic hydrogen. 

Equilibrium fractions of an oxygen beam in various charge states versus 
energy for molecular hydrogen. 



Figure 9 


Figure 10 


Figure 11 


Figure 12 


Figure 13 


Stopping power (or loss function) as a function of oxygen energy for an 
oxygen beam in molecular hydrogen. The nuclear (elastic) contribution is 
shown separately. The stopping power from Northcliffe and Schilling 
(1970) is also shown. 

Relative differential ionization cross sections versus electron energy for 
collisions of oxygen ions at various energies with molecular hydrogen, but 
only for the electron stripping process. 

Energy deposition rate as a function of altitude on Jupiter for incident 
oxygen beams whose energy spectra have the lower energy cut-offs 
indicated. 

Energy deposition rate as a function of altitude for an incident oxygen beam 
whose energy spectrum is a power law (see text) for energies in excess of 
40 keV/amu but with zero flux below this cutoff energy. Contributions of 
different types of collisions (or processes) to the total energy deposition are 
shown. A total atmospheric number density scale(mostly H 2 ) is provided 

on the righthand side of the figure. 

Fractions of the energy deposition contributed by collisions of oxygen in 
different charge states (as labeled) in molecular hydrogen. The total energy 
deposition rate relative to these curves is shown in Figure 12. The sum of 
these fractions for all charge states does not add to 1 because some of the 
energy is being deposited into H and not into H 2 , especially at higher 


altitudes. 



Figure 14 


Figure 15 


Figure 16 


Figure 17 


Figure 18 


Figure 19 


The production rate of secondary electrons versus electron energy, at an 
altitude of 450 km for an incident oxygen spectrum with a 40 keV/amu 
cutoff. 

Production rate of molecular hydrogen ions versus altitude on Jupiter for an 
incident oxygen beam with a 40 keV/amu cutoff. Production from solar 
extreme ultraviolet (EUV) radiation is also included. Contributions to the 
total H 2 + production rate are shown for different processes (primary 

oxygen ions, secondary electrons, and solar EUV radiation). 

Similar to Figure 15 but for the production of H + ions. H + ions can either 
be produced via direct ionization of H by primary oxygen ions, secondary 
electrons, or solar EUV radiation, or by dissociation ionization of H 2 via 

the various types of collisions. 

Differential secondary electron flux directed both up and down as a function 
of altitude for an electron energy of 25 eV. The incident oxygen beam had a 
lower energy cutoff of 40 keV/amu. 

Differential secondary electron flux, at a fixed altitude, as a function of 
electron energy. The electron distribution at this altitude is almost isotropic. 

Number densities of molecular hydrogen, atomic hydrogen, and helium 
versus altitude. The curve labeled unconverged represents H density 
calculated without auroral energy input and only with solar radiation input. 
The curve labeled converged represents the H density calculated with the ion 
aurora turned on, and after convergence of the program. 



Figure 20 Calculated ionospheric profiles (H 3 + and H + ions, and total electron 


density) versus altitude for the ion aurora with "converged" H densities. 
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TABLE 1 


Initial Beam Energies 


k E|<(keV/amu) 


1 

0.3 

25 

230.0 

2 

1.0 

26 

260.0 

3 

2.0 

27 

300.0 

4 

3.0 

28 

340.0 

5 

4.0 

29 

360.0 

6 

5.0 

30 

400.0 

7 

6.0 

31 

440.0 

8 

7.0 

32 

460.0 

9 

8.0 

33 

500.0 

10 

9.0 

34 

530.0 

1 1 

10.0 

35 

560.0 

12 

15.0 

36 

600.0 

13 

20.0 

37 

630.0 

14 

30.0 

38 

660.0 

15 

40.0 

39 

700.0 

16 

50.0 

40 

730.0 

17 

60.0 

41 

760.0 

18 

70.0 

42 

800.0 

19 

80.0 

43 

830.0 

20 

90.0 

44 

860.0 

21 

100.0 

45 

900.0 

22 

130.0 

46 

930.0 

23 

160.0 

47 

960.0 

24 

200.0 

48 

1000.0 



TABLE 2a 


Column Energy Deposition for Unconverged Atmosphere 


Column Production Ratefcm-Ss- 1 ) Energy 


Process 




Energy Input 

9.25x1012 e V 

3.27x1012 eV(35%) 

9.25(1 2)eV 

100.0% 

Backscattered 

Electrons 

— 

5.58(8)eV* 

5.58(8)eV 

» 0.0% 

Lyman Bands 

2.05(10)** 

1.83(10) 

3.88(10) 4.0% 

(65.7kR) 

Werner Bands 

1.51(10) 

1.18(10) 

2.69(10) 

3.0% 

Lyman alpha 

8.97(9) 

5.64(9) 

1.46(10) 

1 .6% 

h 2 + 

1.61(11) 

3.09(10) 

1.92(11) 

33.0% 

H + (H 2 ) 

5.17(10) 

4.37(9) 

5.61(10) 

12.4% 

H + (H) 

5.49(8) 

1.54(8) 

7.03(8) 

0.1% 

Vibrational*** 
Excit. (direct) 

— 

1.57(12) 

1.57(12) 

9.2% 

Vibrational*** 
Excit. (cascade) 

1.59(11) 

1.35(11) 

2.94(11) 

1.7% 

H(H 2 , direct) 

5.17(10) 

1.70(11) 

2.24(11) 

5.1% 

Electron Heat 

— 

2.00(1 1)eV 

2.00(1 1)eV 

2.2% 

Neutral Heat 
(direct) 

3.36(1 1)eV 

6.12(1 1)eV 

9.48(1 1)eV 

10.3% 

01,11,111 Excit. 

9.30(1 1)eV 


9.30(1 1)eV 

10.0% 


Nuclear/Elastic 


(neutral heat) 1.76(1 1)eV 
Misc. - 


1 .76(1 1)eV 1.9% 
5.5% 


100 . 0 % 


* 5.58(8 )eV should be read as 5.58 x 10 8 eV/cm 2 /s 
** 2.05(10) should be read as 2.05 x 10 10 crn^s- 1 
*** Total number of vibrational quanta produced 




TABLE 2b 


Column Energy Deposition for Converged Atmosphere 


Column Production Rate (cm^s* 1 ) 

Energy 

Process Primaries Secondaries Total Efficiency 


Energy Input 

9.25X10 1 * eV 

3.05x1 012©V(33%) 9.25(1 2)eV 100.0% 

Backscattered 

Electrons 

— 

2.08(8) eV* 

2.08(8) eV 

0% 

Lyman Bands 

1.73(10)** 

1.43(10) 

3.16(10) 

3.3% 

Werner Bands 

1.27(10) 

9.43(9) 

(53.7kR) 

2.21(10) 

2.5% 

Lyman alpha 

2.50(10) 

9.21(9) 

3.42(10) 

3.8% 

h 2 + 

1.36(11) 

2.44(10) 

(38.4kR) 

1.60(11) 

27.7% 

H + (H 2 ) 

4.39(10) 

3.49(9) 

4.74(10) 

10.5% 

H + (H) 

2.95(10) 

7.05(9) 

3.66(10) 

8.1% 

Vibrational**’ 
Excit. (direct) 

— 

7.93(11) 

7.93(11) 

4.6% 

Vibrational*** 
Excit. (cascade) 

1.56(11) 

1.02(11) 

2.58(11) 

1.5% 

H(H 2> direct) 

4.39(10) 

1.29(11) 

1.74(11) 

4.3% 

Electron Heat 

— 

5.28(1 1)eV 

5.28(1 1)eV 

5.7% 

Neutral Heat 
(direct) 

2.85(1 1)eV 

4.70(1 1)eV 

7.55(1 1)eV 

8.2% 



01,11,111 Excit. 

9.30(1 1)eV 

— 

9.30(1 1)eV 

10.0% 

Nuclear/Elastic 
(neutral heat) 

1 .70(1 1 )eV 

— 

1.70(1 1)eV 

1.8% 

Misc. 

— 

— 


8.0% 

100.0% 


* 2.08(8)eV should be read as 2.08 x 10 8 eV/cm 2 /s 
1.73(10) should be read as 1.73 x 10 10 cm^s* 1 
Total number of vibrational quanta produced 




TABLE 3 


Comparison of Energy Efficiencies for Ion and Electron 
Precipitation into the Unconverged Atmosphere for Selected 
Processes 

Process Oxygen (E>4 0 keV/amu l 10 keV electrons 


H 2 + 

33.0% 

38.9% 

H+(total) 

12.5% 

2.8% 

Lyman/Werner bands 

7.0% 

14.5% 

01,11,111 excitation 

10.0% 

0.0% 

(all vibration + direct 



neutral heat + electron 



heat) 

23.4% 

24.4% 

direct H production 

5.1% 

6.1% 


